Dynamical deformation effects in subbarrier fusion of ^^Ni+^^^Sn 
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We show that dynamical deformation effects play an important role in fusion reactions involving 
the ^*Ni nucleus, in particular the ®''Ni+^''^Sn system. We calculate fully microscopic interaction 
potentials and the corresponding subbarrier fusion cross sections. 

PACS numbers: 21.60.-n,21.60. Jz 
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Recently observed enhanced subbarrier fusion cross 
sections for the neutron rich ^''Ni+^'^^Sn system has 
further invigorated the research in low-energy nuclear 
reactions involving exotic nuclei. For this system fu- 
sion cross-sections were measured in the energy range 
142 MeV < Sc.m. < 176 MeV. In particular, it was 
found that fission is negligible for i?c.m. < 160 MeV 
and therefore the evaporation residue cross-sections have 
been taken as fusion cross sections. The enhancement of 
subbarrier fusion was deduced from comparison with a 
barrier penetration calculation, using a phenomenologi- 
cal Woods-Saxon interaction potential whose parameters 
were fitted to reproduce the evaporation residue cross sec- 
tions for the ^^Ni-l-^^^'Sn system. Similarly, sophisticated 
coupled-channel calculations 0, Q , which are known to 
enhance the fusion cross sections by considering coupling 
to various excitation channels and neutron transfer, have 
significantly underestimated the subbarrier fusion cross 
sections for the ^"^Ni-f -'^^^Sn system 1]. 

In general, the fusion cross sections depend on the in- 
teraction potential and form factors in the vicinity of the 
Coulomb barrier. These are expected to be modified dur- 
ing the collision due to dynamical effects. In addition, 
experiments on subbarrier fusion have demonstrated a 
strong dependence of the total fusion cross section on 
nuclear deformation 0|. The dependence on nuclear ori- 
entation has received particular attention for the forma- 
tion of heavy and superheavy elements Q and various 
entrance channel models have been developed to predict 
its role in enhancing or diminishing the probability for 
fusion HQ. 

Recently, we have developed a new approach for cal- 
culating heavy-ion interaction potentials which incorpo- 
rates all of the dynamical entrance channel effects in- 
cluded in the time-dependent Hartree-Fock (TDHF) de- 
scription of the collision process Q . These effects include 
the neck formation, particle exchange, internal excita- 
tions, and deformation effects to all order, as well as the 
effect of nuclear alignment for deformed systems 0, ■ 
The method is based on the TDHF evolution of the nu- 
clear system coupled with density-constrained Hartree- 
Fock (DCHF) calculations % J^] to obtain the interac- 
tion potential, given by 



(1) 



tion, such as the Skyrnie force, V^K) can be constructed 
by performing a TDHF evolution and minimizing the en- 
ergy at certain times to obtain Edc{R), while the nuclear 
binding energies Ea-^ and Ea2 are the results of a static 
Hartree-Fock calculation with the same effective interac- 
tion [|. 

We have carried out a number of TDHF calculations 
with accompanying density constraint calculations to 
compute V(R) given by Eq. A detailed description 
of our new three-dimensional unrestricted TDHF code 
has recently been published in Ref flU . For the effective 
interaction we have used the Skyrme SLy5 force [T^ in- 
cluding all of the time-odd terms. In our case the ^'^Ni 
nucleus is essentially oblate with a small mix of triaxial- 
ity, having a quadrupole moment of -0.45 b. This is also 
confirmed by other calculations ^3 Q| and suggested by 
experiments TDHF calculations were initialized at 
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FIG. 1: (Color online) Internuclear potential obtained from 
Eq. Q for the head-on collision of the ^''Ni-l-^^^Sn system 
at -Ec.m. = 158 MeV. The dashed line shows the empirical 
Woods-Saxon potential used in Ref. 
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158 MeV. In Ref. 



we have shown that the 



The potential deduced from Eq. contains no param- 
eters and no normalization. Given an effective interac- 



calculation of the potential barrier is not sensitive to the 
choice of the TDHF initialization energy, the only differ- 
ence being a slightly lower potential well for lower ener- 
gies. In Fig. n we show the results obtained for the inter- 
action potential as well as the empirical potential barrier 
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used in Ref. [jj]. The angle (3 indicates the orientation of 
the symmetry axis. In the case oi (3 — Q° the symmetry 
axis of the oblate ^^Ni is aligned with the collision axis 
and for (3 = 90° the symmetry axis is perpendicular to 
the collision axis. For the case of parallel orientation the 
calculated barrier is almost exactly the same as the one 
used in Ref. 0, having a barrier height of 155.8 MeV. 
The difference for smaller R values is due to the use of 
the point Coulomb interaction in the model calculation, 
which is unphysical when nuclei overlap. The same ar- 
gument applies to small differences at large R values, 
since the Coulomb interaction is slightly different due to 
the deformed Ni nucleus. We would like to emphasize 
again that our calculations do not contain any param- 
eters or normalization. On the other hand, the barrier 
corresponding to the perpendicular alignment is consid- 
erably lower, peaking at 150.1 MeV, and has a narrower 
width. 

The physical picture which emerges from these calcu- 
lations is that for center-of-mass energies in the range 
150.1-155.8 MeV the fusion cross section would be dom- 
inated by the channel above the lower barrier since the 
contribution via tunneling through the higher barrier will 
be substantially smaller. Similarly, for energies below 
150.1 MeV transmission through the lower barrier will 
produce the dominant contribution. Of course, for ener- 
gies above 155.8 MeV both barriers will contribute. As 
a result, the only data point which is truly subbarrier is 
the lowest energy point. In Fig.|21we show the calculated 
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FIG. 2: (Color online) Fusion cross sections obtained for 
the "^^Ni-f^^Sn system using the microscopically calculated 
potentials discussed in the manuscript. Also shown (circles) 
are the experimental values from Ref. Q. 



fusion cross sections as a function of the center-of-mass 
energy. We have used a simple WKB approach to calcu- 
late the cross section for the lowest energy point and the 
parabolic approximation via the Wong formula for 
the higher energy points. Also shown are the correspond- 
ing experimental values (circles) as well as the barrier 
peneration model results (dashed line) from Ref. P]. As 
anticipated, the calculated cross sections, with the excep- 
tion of the lowest energy data point, are now above the 
corresponding experimental values. There are a number 
of reasons for over-predicting the data at higher energies. 
The first reason is the alignment probability of the de- 
formed nucleus, which could be calculated with the avail- 
ability of the excitation spectrum for the ^^Ni nucleus 0] . 
Physically, one expects a distribution of barriers starting 
from the lowest barrier and approaching the highest bar- 
rier. The second factor is the quality of the parabolic ap- 
proximation used in the Wong formula. It is well known 
that the rising Coulomb tail of barriers cannot be prop- 
erly accounted for by a single parabola, thus resulting in 
a somewhat thinner barrier and a larger fusion cross sec- 
tion. In addition, for energies above 160 MeV the fission 
channel opens up. Finally, despite improving the lowest 
energy cross section by many orders of magnitude we find 
a cross section of 0.0035 mb, which is still a factor of 200 
lower than the experimental value. 



In conclusion, we have performed microscopic calcu- 
lations of the interaction potentials for the ^^Ni-l-^'^^Sn 
system. We observe that dynamical deformation effects 
play a very significant role in the calculation of fusion 
cross sections. This observation further underscores the 
necessity of detailed structure information for neutron 
and proton rich systems for the better description of fu- 
sion cross sections involving these nuclei. Availability 
of detailed structure data for the ^^Ni nucleus may help 
explain the discrepancy for the lowest energy point. 



This work has been supported by the U.S. Depart- 
ment of Energy under grant No. DE-FG02-96ER40963 
with Vanderbilt University. Some of the numerical cal- 
culations were carried out at the IBM-RS/6000 SP su- 
percomputer of the National Energy Research Scientific 
Computing Center which is supported by the Office of 
Science of the U.S. Department of Energy. 



[1] J. F. Liang et al, Phys. Rev. Lett. 91, 152701 (2003); [2] A. B. Balantekin and N. Takigawa, Rev. Mod. Phys. 70, 
Phys. Rev. Lett. 96, 029903(E) (2006). 77 (1998). 



3 



[3] H. Esbcnscn, Prog. Thcor. Phys. Suppl. 154, 11 (2004). 
[4] R. G. Stokstad, Y. Eiscn, S. Kaplanis, D. Pcltc, U. Smi- 

lansky, and I. Tscrmya, Phys. Rev. Lett. 41, 465 (1978). 
[5] K. Nishio, e£ al, J. Nucl. Radiochemical Sci., 3, 89 

(2002). 

[6] C. Simenel, Ph. Chomaz, and G. de Prance, Phys. Rev. 

Lett. 93, 102701-1 (2004). 
[7] A. S. Umar and V. E. Oberacker, Phys. Rev. C 74, 

024606 (2006). 

[8] A. S. Ulnar and V. E. Oberacker, Phys. Rev. C 74, 
021601(R) (2006). 

[9] R. Y. Cusson, P. -G. Reinhard, M. R. Strayer, J. A. 
Maruhn, and W. Greiner, Z. Phys. A 320, 475 (1985). 
[10] A. S. Umar, M. R. Strayer, R. Y. Cusson, P. -G. Rein- 
hard, and D. A. Bromley, Phys. Rev. C 32, 172 (1985). 



[11] A. S. Umar and V. E. Oberacker, Phys. Rev. C 73, 

054607 (2006). 

[12] E. Chabanat, P. Bonchc, P. Hacnscl, J. Meyer and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998); Nucl. Phys. 

A643, 441 (1998). 
[13] J. Dobaczewski, M. V. Stoitsov, and W. Naaarewicz, 

Skyrme-HFB deformed nuclear mass table, ed. R. Bijker, 

R.F. Gaston, and A. Prank, (AIP, New York, 2004), 726, 

51 (2004). 

[14] G. A. Lalazissis, S. Raman, and P. Ring, At. Data Nucl. 

Data Tables 71, 1 (1999). 
[15] J. J. Vega, E. F. Aguilera, G. Murillo, J. J. Kolata, A. 

Morsad, and X. J. Kong, Phys. Rev. C 42, 947 (1990). 
[16] G. Y. Wong, Phys. Rev. Lett. 31, 766 (1973). 



